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In response to recent nuclear-magnetic-resonance measurements on a ferrimagnetic chain com- 
pound NiCu(C7H 6 N206)(H 2 0)3-2H20 [Solid State Commun. 113 (2000) 433], we calculate the 
nuclear spin- lattice relaxation rate 1/Ti in terms of a modified spin- wave theory. Emphasizing 
that the dominant relaxation mechanism arises from the direct (single-magnon) process rather than 
the Raman (two-magnon) one, we explain the observed temperature and applied-field dependences 
of 1/Ti. Ferrimagnetic relaxation phenomena are generally discussed and novel ferrimagnets with 
extremely slow dynamics are predicted. 

PACS numbers: 75.10.Jm, 75.30. Ds, 76.60.Es 



I. INTRODUCTION 

Design of molecule-based ferromagnets has been one of 
the most exciting subject in Materials science. One can 
in principle obtain molecular ferromagnets by assembling 
molecular bricks so as to construct a low-dimensional sys- 
tem with a magnetic ground state and then coupling the 
chains or the layers again in a ferromagnetic fashion jj] . 
In the naivest attempt to obtain a magnetic ground state, 
we may couple the nearest-neighbor magnetic centers fer- 
romagnetically g|. However, it is often difficult to realize 
the symmetry conditions favoring the parallel alignment 
of local spins. The difficulty was overcome by the intro- 
duction of a new concept— antiferromagnetically coupled 
polymetallic systems with irregular spin-state structures 
H . Ordered bimetallic chain compounds were thus syn- 
thesized and since then the magnetic properties of ferri- 
magnetic Heisenberg chains have extensively been inves- 
tigated @-|io|. 

Recently nuclear magnetic resonance for 1 H nuclei has 
been performed in a ferrimagnetic chain compound and 
stimulative temperature and applied-field dependences of 
the proton spin-lattice relaxation rate 1 jT\ has been re- 
ported pl| . The authors analyzed 1/Ti in terms of the 
naivest spin-wave theory jl2| but could not successfully 
interpret the characteristic field dependence, which looks 
like 1/Ti cx l/\/~H. We here point out that their ar- 
gument broke down because they attributed the dominant 
relaxation mechanism to the Raman (two-magnon) pro- 
cess. Modifying the naivest spin-wave theory so as to 
fully describe the thermodynamics, we demonstrate that 
the direct (single-magnon) relaxation process can be ef- 
fective in Heisenberg ferrimagnets. We further discuss 
potential ferrimagnetic spin dynamics arising from the 
characteristic twofold excitations |lj|. 



II. SPIN- WAVE APPROACH 

The measured compound NiCu(pba)(H 2 0)3-2H 2 
(pba = 1, 3-propylenebis(oxamato) = C7H 6 N 2 06) JTJ] 



consists of ordered bimetallic chains with alternating oc- 
tahedral Ni 2+ and square-pyramidal Cu 2+ ions bridged 
by oxamato groups. The one-dimensional character holds 
down to 7[K] under the exchange coupling J//cb — 
121 [K]. The g factors of the S = 1 and s — \ spins 
are both close to 2 |15|j . Thus the material is reasonably 
described by the one-dimensional mixed-spin Heisenberg 
Hamiltonian 



N 
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H = JJ2 (Si ■ s 3 + Sj ■ S j+1 ) - 9m Hj2(S- +8$, 

(1) 

where Sj and Sj are respectively spin-S'(= 1) and spin- 
s' operators at the jth elementary cell. The spin- 
lattice relaxation rate is given by 
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where wn = 7N-ff is the Larmor frequency of the nuclei 
with the gyromagnetic ratio 7n and the summation ^ 
is taken over all the electronic eigenstates \n) with energy 
E n . Assuming the hyperfine interaction to be isotropic, 
we may represent Hj n t as 



H m , = Y,{AI r S 3 + BI r s 



3) ' 



(3) 



with the nuclear spin operators Ij and the dipolar cou- 
pling constants A and B. 

In order to calculate 1/Ti in practice, we introduce the 
bosonic operators for the spin deviation as 
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at = b] {2s - b]bj)^ 2 , s) = -s + b]bj 



(4) 



Then, assuming O(S) — 0(s), we can expand Hint with 
respect to 1/S and obtain 1/Ti = J2i=i 2 ••■ V^i 
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where 1/T X W is the I -magnon relaxation rate within the 
first-order spin-lattice relaxation process and is given by 
replacing H int by its 0(S' 1_ '/ 2 )-component in Eq. (||). 
On the other hand, within the up-to-0(S ) approxima- 
tion, the electronic Hamiltonian (|l|) can be diagonalized 
as M 



H * H SW = E e + fe4o* + 



(5) 



Here, Eg is the spin-wave ground-state energy, whereas 
a\ and /?| create the spin waves of ferromagnetic and 
antiferromagnetic aspects Jl8| , respectively, and are re- 
lated with the sublattice bosons via at = a^cosh^ + 



bkS'mhOk and fr k = afcSinMfe + fcj.coshflfc with a\ = 
(l/ViV^.e-^-V^at, b\ = (l/^Ej^'+^fc}, 

and tanh(26» fe ) = 2v / £scos(fe/2)/(S + s). We here take 
twice the lattice constant as unity. The dispersion rela- 
tions are given by 
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FIG. 1. Dispersion relations of the linear- (dotted lines; 
LSW) and interacting- (solid lines; ISW) spin-wave excita- 
tion energies compared with the lowest-energy bands in the 
subspaces of X^i^ 2 "h s ?) ~ ^ S )N T 1 estimated by a 
quantum Monte Carlo method (x; QMC) for (S, s) = (1, |). 



The linear spin waves are characterized by ajf , while the 
interacting spin waves by ui^ with the O(S ) quantum 
corrections Sui^. We plot ui^ and uf in Fig. |l]. The 
spin waves describe the elementary excitations satisfac- 
torily, especially the low-energy excitations near the zone 
center, which must be the most relevant to the nuclear 
spin relaxation. The antiferromagnetic excitation ener- 
gies are underestimated by the linear spin waves but are 
significantly corrected in consideration of the interactions 
between them. All the following calculations are based 
on the interacting-spin-wave dispersions ui^ . 

The 5 function in Eq. (j^) insures the conservation of 
energy in the transition. Considering the significant dif- 
ference between the electronic and nuclear energy scales 
(Tilj-n 10~ 5 J), the relevant spin-wave excitations are 
strongly limited. Only the small-momentum ferromag- 
netic spin waves contribute to the direct process. In 
the Raman process, both ferromagnetic and antiferro- 
magnetic spin waves are effective but interband transi- 
tions are still irrelevant. Scatterings between ferromag- 
netic and antiferromagnetic spin waves are relevant to the 
higher-order processes containing three or more magnons. 
We show a few leading terms of the relaxation rate: 



4( 5 MB?i7N) 2 / 2(S - s) 



r« njh y Sshu N /j 

x (y/SAcosh9 - y/sBsmhe ) 2 {riQ + 1) 
1 ^ 4( 5A1b ?>7n) 2 (S - s) 



(12) 



T (2) ~ NJh^iSsk) 2 + 2{S - s)SshLU N /J 

x [(Acosh 2 ^ - Ssinh 2 ^)^^ + 1) 

+ (Asinh 2 e k - Bcosh 2 6 k ) 2 n+{n+ + 1)] , 

_1 \p (9^b^1n) 2 (S - s) 

^(3) ~ 2^ 



E 



(13) 



2N 2 Jh^jS 2 s 2 {k 2 + k' 2 ) + 2(5 - s)Sshu N /J 



[( 



A 



x | ( ^^cosh^cosMfe/CosMfcifc/ 
'V5 



D 



=sinh(9fc 



+ 4 



. A 



where n. 



(a{a k ) 



1 2 _ 



and ri£ 



k>( n k+k> + 


1) 






h+k' 7= 


sinh^ 






l k>( n k+k> + 


1)" , 



(14) 



(plPk) are the distri- 



bution functions of the spin waves. In general 1/T-f 1 
rapidly decreases as a function of I. However, a slight 
anisotropy, for example, turns off the direct process and 
therefore makes the Raman process effective. In such a 
case, the three-magnon-process relaxation rate can also 
be effecitve, being enhanced via the exchange interaction, 
at high temperatures ]16|,|l9|| . 
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Now the problem is reduced to the evaluation of rvt. 
The naivest thermodynamics defining the partition func- 
tion as Z = Tr[e~ Wsw / fcBT ] ends in the divergence of 
with increasing temperature. Thus we here employ the 
modified spin wave theory for ferrimagnets ]r^ , |20| and 
obtain the spin-wave distribution functions as 
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where /i is a Lagrange multiplier controlling the staggered 
magnetization and is self-consistently determined under 
the constraint ^2 k J2 a =± n h/ U} 'k = N. We show in Fig. ^ 
the thus-obtained as functions of k. At low tempera- 
tures, exhibits a pronounced peak at k = 0, whereas 
n\ stays negligibly small in comparison with , because 
the ferromagnetic and antiferromagnetic excitations are 
respectively gapless and gapped from the ground state. 
Even at high temperatures, remains less dominant 



due to the large gap w^L = 2J(S— s)(l + r 2 /\/Ss) with 
r 2 = 0.4777 for (S, s) = (1, §). The applied field, slightly 
shifting the spin- wave excitation energies, reduces 
and enhances n^. 
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FIG. 2. The momentum distribution functions of the ferro- 
magnetic (a) and antiferromagnetic (b) spin waves at various 
values of temperature and the applied magnetic field in the 
case of (S, s) = (1, |). 



ture to the relaxation rate. Therefore, l/X^ results in a 

(2) 

monotonically decreasing function of T, whereas l/T^ 
contains increasing components as well as overwhelming 
decreasing ones. Here is a fascinating parameter to be 
adjusted, r = A/ B. Dipolar coupling constants are pro- 
portional to the inverse cube of the distance between the 
interacting nuclear and electronic spins and therefore r 
is quite sensitive to the crystalline structure. Consid- 
ering the predominance of the small-momentum excita- 
tions, we find there is a special point of r = s/S, where 
the dominant relaxation mechanism mediated by the fer- 
romagnetic spin waves does not work at all in every l- 
magnon process and thus the Raman process mediated 
by the antiferromagnetic spin waves becomes the lead- 
ing relaxation mechanism. That is why as r approaches 
s/S, the Raman contribution increases emphasizing the 
increasing behavior as a function of T . 



7500 



6000 



4500 



3000 



1500 




3. 



200 



150 \ 



100 



50 



(a) r = 0.0 Direct Direct+Raman 

0.5 [T] 

1.0 [T] 
3.0 [T] 
5.0 [T] 



S (b) r = 0.4 

\ 
\ 
\ 

!\ \ 

: \ x 



1 2 
kT/J 



FIG. 3. Temperature dependences of the direct-process 
and Raman-process relaxation rates at various values of the 
applied magnetic field in the case of (S, s) — (1, |). l/T^ is 



plotted by thin lines, while + 1/T^> by thick lines 



-,(2) 



III. EXPLANATION OF THE EXPERIMENT 



Let us observe general behavior of 1/Ti in Fig. ||. We 
have learned in Fig. |2| that the ferromagnetic and anti- 
ferromagnetic spin waves respectively contribute decreas- 
ing and increasing components as functions of tempera- 



We apply the theory to the proton spin relaxation 
jll| in NiCu(pba)(H 2 0)3-2H 2 0. Since the protons mainly 
contributing to 1 jT\ turn out to lie in the pba groups near 
Cu ions, we hear set r smaller than unity. 
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FIG. 4. Temperature dependences of the direct-process re- 
laxation rate plus the Raman-process one, + 1/T^ 2 \ 
at various values of the applied magnetic field and the struc- 
tural parameters in the case of (S, s) — (1, ~), assuming that 
J/k-B == 121 [K]. The calculations are shown by lines, where 
H = 0.53, 1.63, 3.15[T] from the top to the bottom, while the 
measurements by symbols. 
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FIG. 5. Field dependences of the direct-process relax- 



ation rate plus the Raman-process one, 1/t[ + , at 

various values of the structural parameters in the case of 
(S,s) — (1, |), assuming that J//cb = 121[K]. Temperature 
is set equal to 280 [K]. The calculations are shown by lines, 
while the measurements by o. A calculation for 1/t[ 
shown for reference. 
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is also 



We analyze the observed temperature dependences of 
1/Ti [[llj in Fig. |[ Though no parametrization ends 
in a total agreement with the measurements, the calcu- 
lation can reproduce them fairly well. In the present 
sample, the relevant protons have a rather wide distribu- 
tion around Cu ions and therefore the distances between 
the interacting proton and electron spins, which we refer 
to as a and b for S — 1 and s — i, respectively, can not 
be determined definitely. If we assume that A ~ 1/a 3 

and B ~ 1/b 3 , 4.48 x 1(T 3 [A _ ] for B may be too small, 
while 0.4 for r be too large. As the applied field increases, 



the recovery curve of the spin-echo decay more and more 
deviates from a single-exponential function pjj . The ex- 
perimental estimates of 1 jT\ at H > 1 [T] seem to contain 
larger uncertainty. 

Field de pen dence of 1 /Xi was measured in more detail 
at 280[K] |l|] and is analyzed in Fig. |. We again find 
that the calculations with small r can generally be fitted 
to the measurements. We emphasize that the linear de- 
pendence of 1/Ti on \j\fH with a steep slope can never 
be attributed to the Raman process but be described 
in consideration of the direct process. As the most sig- 
nificant field dependence is obtained at r = in each 
relaxation process, the Raman contribution to the field 
dependence turns out much smaller than that from the 
direct process. The three-magnon contribution, even if 
it is enhanced via the exchange interaction, is still less 
relevant. In general, the momentum integral in the cal- 
culation of 1/Ti weakens the linear field dependence and 
even leads to logarithmic one at high fields. 

We should be reminded that the present spin-wave the- 
ory, which is free from both quantum and thermal diver- 
gences of the sublattice magnetizations, is highly suc- 
cessful but its ability to describe thermodynamics is still 
restricted quantitatively flS| ]. Considering minor correc- 
tions to the theory as well, arising from distribution of 
proton spins, slight difference between g factors, weak 
interchain interactions and so forth, we are allowed to 
conclude that the direct process plays a leading role in 
the nuclear spin relaxation of NiCu(pba)(H20)3-2H20. 



IV. CONCLUDING REMARKS 

The dominant relaxation mechanism for the title com- 
pound has been shown to arise from the direct process. It 
is at the same time the demonstration that the material 
is a really isotropic magnet. 

Another harvest we have obtained is the prediction of 
novel ferrimagnets with extremely slow dynamics. Un- 
der the crystalline structure with r = s/S, the fast re- 
laxation processes mediated purely by the ferromagnetic 
spin waves hardly work and thus the nuclear spin is very 
slow in relaxing especially at low temperatures. 
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